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Abstract In this Chapter, we present recent theoretical developments on the fi- 
nite temperature transport of one dimensional electronic and magnetic 
quantum systems as described by a variety of prototype models. In 
particular, we discuss the unconventional transport and dynamic - spin, 
electrical, thermal - properties implied by the integrability of models as 
the spin- 1/2 Heisenberg chain or Hubbard. Furthermore, we address 
the implication of these developments to experimental studies and the- 
oretical descriptions by low energy effective theories. 
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1. Introduction 

The electronic and magnetic properties of reduced dimensionality ma- 
terials are significantly modified by strong correlation effects. In par- 
ticular, over the last few years, the physics of quasi-one dimensional 
electronic systems, has been the focus of an ever increasing number of 
theoretical and experimental studies. They are realized as three dimen- 
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sional (3D) compounds composed of weakly interacting chains or, in 
a latest and very promising development, as monoatomic width chains 
fabricated by self-assembly on surfaces. 

Experimentally, recent studies made possible by the synthesis of new 
families of compounds characterized by very weak interchain coupling 
and low disorder, indicate unconventional transport and dynamic be- 
havior; for example, unusually high thermal conductivity in quasi-one 
dimensional magnetic compounds [1, 2, 3], ballistic spin transport in 
magnetic chains [4, 5] or optical conductivity in quasi-lD organic con- 
ductors, showing a low frequency very narrow "Drude peak" even at 
relatively high temperatures [6, 7]. 

Theoretically, it is well known that one dimensional (ID) systems of 
interacting electrons do not follow the phcnomcnological description of 
the ordinary Landau Fermi liquids, but rather they are characterized 
by a novel class of collective quantum states coined Luttinger liquids 
[8]. Furthermore, it has quite recently been realized that several proto- 
type models commonly used to describe ID materials imply ideal trans- 
port properties (dissipationless) even at high temperatures. This phe- 
nomenon is the quantum analogue of transport by nondecaying pulses 
(solitons) in ID classical nonlinear integrable systems [9]. 

Of course, ID electronic and magnetic systems have, since the sixties, 
been a favorite playground where theoretical ideas were confronted with 
experimental results on an ever improving quality and variety of quasi- 
lD compounds. We are now in a position to claim reliable theoretical 
analysis on the thermodynamics, quantum phase transitions and spectral 
functions of prototype many body Hamiltonians used to describe ID 
materials. The tools at hand range from exact analytical solutions (e.g. 
using the Bcthc ansatz (BA) method) [10], the low energy Luttinger 
liquid approach and powerful numerical simulation techniques as the 
Quantum Monte Carlo (QMC) [11] and Density Matrix Renormalization 
Group (DMRG) [12] method. In particular, the ground-state properties 
as well as the low-tcmpcraturc behavior of the correlation functions and 
of most static quantities in the scaling (universal) regime of Luttinger 
liquids have been extensively studied in recent decades and are well 
understood by now [13, 14, 15, 16]. 

On the other hand, although most experimentally relevant, less stud- 
ied and understood are the transport and dynamic properties of ID 
interacting electronic or magnetic systems. Paradoxically, while the 
equilibrium properties of prototype integrable models as the spin-1/2 
Heisenberg or Hubbard model have been extensively analyzed, their fi- 
nite temperature transport has attracted little attention; most studies 
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till a few years ago, have basically focused on the low energy descrip- 
tion in terms of the Luttinger model. A noticeable exception has been 
the issue of diffusive versus ballistic behavior and thermal conductivity 
of spin chains, a long standing and controversial issue. The difficulties 
encountered with the transport quantities can be attributed to the fact 
that the scattering and dissipation in clean ID fermionic systems are not 
dominated by low-energy processes and thus the transport properties are 
not universal. 

Presently, transport properties are at the focus of intense theoretical 
activity; in particular, prototype integrable models (as the Heisenberg, t- 
J, Hubbard, nonlinear-o") are studied by exact analytical techniques (e.g. 
Bethe ansatz, form factor method) and numerical simulations. However, 
the complexity of these methods often renders the resulting behavior 
still controversial. Furthermore, the transport of quasi-one dimensional 
systems is (re-) analyzed within the effective Luttinger liquid theory or 
by semiclassical, Boltzmann type, approaches. 

It is fair to say that the study of finite temperature/frequency con- 
ductivities in strongly correlated systems presents at the moment funda- 
mental conceptual as well as technical challenges. Development of new 
analytical and numerical simulation techniques is required, as well as 
progress on the basic understanding of scattering mechanisms and their 
effects. 

In the following, we will mostly concentrate on the conductivity of 
bulk, clean systems where the scattering mechanism is due to electronic 
or magnetic interactions (Umklapp processes). In particular, we will not 
address issues on the transport of mesoscopic systems (e.g. nanowires, 
nanotubes) or other dissipation mechanisms as coupling to phonons or 
disorder. 

In section 2 we start by presenting some elements of linear response 
theory (or Kubo formalism), the theoretical framework commonly used 
for describing transport properties. Then, in section 3 and 4 wc continue 
with a presentation of the state of the transport properties of prototype 
systems, in particular the Heisenberg and Hubbard model. In section 
5, we present a short overview of alternative approaches based on low 
energy effective field theories as the Luttinger liquid, sine-Gordon and 
quantum nonlinear-cr models. Finally, in section 6, we close with a 
critical assessment of the present status and a discussion of open issues. 

This presentation is definitely not an exhaustive account of theoretical 
studies on the transport properties of one dimensional quantum systems 
but it rather aims at presenting a coherent and self-contained view of 
some recent developments. 
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2. Linear response theory 

In this section we introduce the basic definitions and concepts that 
will be used in the later development. The framework of most trans- 
port studies is linear response theory where the conductivities arc given 
in terms of finite temperature (T) dynamic correlations calculated at 
thermodynamic equilibrium [17]. For instance, the real part of the elec- 
trical conductivity at frequency lo is given by the corresponding dynamic 
current correlation Xjj(a;), 



a'{L0) = 2-nD5{uj) + (Treg{oj) (1.1) 
Gregiu^) = ^J-Xjji^^) (1-2) 

XiiM = i dte'''{[j{t),j]), z = u + irj (1.3) 
Jo 

with j the appropriate current operator. In a spectral representation 
the conductivity is, 

1 — e"^'^ TT 

CTregii^) = y V^Pn Y] I I j I ™) P ^(^ " (^m - En)), (1-4) 

I n),en denoting the eigenstates and eigenvalues of the Hamiltonian, p„ 
the corresponding Boltzmann weights and /3 the inverse temperature (in 
the following we take h = kb = e = 1); the dc conductivity is given by 
the limit adc = c^regi^ 0). We will mostly discuss one dimensional 
tight-binding models on L sites where the current operator does not 
commute with the Hamiltonian. 

To define the current operators we use the continuity equations of 
charge, magnetization or energy for the electrical, magnetic and thermal 
conductivity, respectively. We will explicitly present them below in the 
discussion of the Heisenberg and Hubbard models. 

A quantity that presently attracts particular attention is the prefac- 
tor D of the 5— function, named the Drude weight or charge stiffness. 
This quantity was introduced by W. Kohn in 1964 as a criterion of 
(ideal) conducting or insulating behavior [18] at T = in the context 
of the Mott-Hubbard transition. This meaning becomes clear by noting 
that D is also the prefactor of the low frequency, imaginary (reactive - 
nondissipative) part of the conductivity, 



D = ^[u;a"(a;)]_o = \ \{-T) - Y.Pn E ' ' ' ' '"^ ' ] ; (1-5) 



5 



here (T) denotes the thermal expectation value of the kinetic energy, 
generalizing the zero temperature expression to T > by consider- 
ing a thermal average. Thus, a finite Drude weight implies an "ideal 
conductor" , a freely accelerating system. The second definition of the 
Drude weight follows from the familiar optical sum-rule [19, 20, 21] using 



with the average value of the kinetic energy replacing for nearest neigh- 
bor hopping tight binding models the usual ratio of density over mass 
of the carriers for systems in the continuum. 

At r = the Drude weight Dq = D{T = 0) is the central quantity 
determining charge transport. As already formulated by Kohn in a very 
physical way, Dq can also be expressed directly as the sensitivity of the 
ground state energy eo to an applied flux (f) = eA (e = 1), 



For a clean system, since at T = there cannot be any dissipation, one 
expects that (Jreg{^ — > 0) = and we have to deal with two fundamen- 
tally different possibilities with respect to Dq: 

-Do > is characteristic of a conductor or metal, 
Dq = characterizes an insulator. 

The insulating state can originate from a filled electron band (usual band 
insulator) or for a non-filled band from electron correlations, that is the 
Mott-Hubbard mechanism; the latter situation is of interest here. Note, 
that the same criterion of the sensitivity to flux has been applied to 
disordered systems, in the context of electron localization theory [22]. 

The theory of the metal-insulator transition solely due to Coulomb 
repulsion (Mott transition) has been intensively investigated in the last 
decades by analytical and numerical studies [23] of particular models of 
correlated electrons and it is one of the better understood parts of the 
physics of strongly correlated electrons. 

At finite temperatures, within the usual Boltzmann theory for weak 
electron scattering, the relaxation time approximation represents well 
the low frequency behavior. 



eq.(l.l). 




(1.6) 



Dq = 



1 dho 
2L '-^-^O" 



(1.7) 



a{uj) = (7dc/(l + ^wr), 



(1.8) 
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where the relaxation time r depends on the particular scattering mech- 
anism and is in general temperature dependent. In the following, we 
consider only homogeneous systems without any disorder, so the rele- 
vant processes in the solid state are electron-phonon scattering and the 
electron-electron (Coulomb) repulsion. When the latter becomes strong 
it is expected to dominate also the transport quantities. 

Even in a metal with Dq > it is not evident which is the relevant 
scattering process determining t{T) and addT). In the absence of disor- 
der and neglecting the electron-phonon coupling the standard theory of 
purely electron-electron scattering would state that one needs Umklapp 
scattering processes to obtain a finite r. That is, the relevant electron 
Hamiltonian includes the kinetic energy H^in, the lattice periodic po- 
tential V and the electron-electron interaction Hint, 

H = Hk^n + V + Hint. (1.9) 

Then, in general, the electronic current density j is not conserved in 
an Umklapp scattering process as the sum of ingoing electron momenta 
equals the sum of outgoing ones only up to a nonzero reciprocal vector 
G, J2i ki = mG. In other words, the noncommutativity of the current 
with the Hamiltonian, [H,j] ^ 0, leads to current relaxation and thus, 
by the fluctuation-dissipation theorem, to dissipation. The interplay 
of V and Hint, however, turns out to be fairly involved in the case of 
strong electron-electron repulsion. This will become clear in examples 
of integrable tight binding models of interacting systems that we will 
discuss below, which have anomalous (diverging) transport coefficients. 

Experiments on many novel materials, - strange metals - with corre- 
lated electrons, question the validity of the concept of a current relax- 
ation rate 1/r. Prominent examples are the superconducting cuprates 
with very anisotropic, nearly planar, transport [23] where the experi- 
mentally observed (t{u)) in the normal state can be phenomenologically 
described only by strongly frequency (and temperature T) dependent 
t{uj, T). The experiments on (t{uj) in quasi-lD systems are covered else- 
where [7]. 

With this background, we will now discuss different possible scenaria 

for the behavior of the T > conductivity. A clean metallic system at 
T = is characterized by a 5— function Drude peak and a finite frequency 
part that vanishes, typically with a power law dependence, implying zero 
dc regular conductivity. In the common sense scenario, at finite temper- 
atures the (5— function broadens to a "Drude peak" of width inversely 
proportional to a characteristic scattering time and thus a finite a; — > 
limit implying a finite dc conductivity. The scattering mechanisms can 




Figure 1.1. Schematic representation 
of the typical behavior of the conduc- 
tivity for a clean metal. 



Figure 1.2. Illustration of the con- 
ductivity of a clean metal remaining 
an ideal conductor at finite tempera- 
tures. 



be intrinsic, due to interactions, or extrinsic due to coupling to other 
excitations, phonons, magnons etc. This typical behavior is shown in 
Fig. 1.1. 

Actually, for a finite size system (as often studied in numerical sim- 
ulations) D is nonzero even at finite T; it only goes to zero, typically 
exponentially fast, as the system size tends to infinity. Physically, this 
expresses the situation where the thermal scattering length is less than 
the system size. 

But it is also possible that constraints on the scattering mechanisms 
limit the current decay, so that the system remains an ideal conduc- 
tor {D > 0) even at T > 0. A schematic representation of a system 
remaining an "ideal conductor" at finite T is shown in Fig. 1.2. 

In a system with disorder, D vanishes even at zero temperature and 
the dc residual conductivity is finite (provided the disorder is not strong 
enough to produce localization). 

For an insulating system, e.g. due to interactions or the band struc- 
ture as we discussed above, D vanishes at zero temperature; in the 
conventional case, D remains zero at T > 0, while activated carriers 
scattered via different processes give rise to a finite dc conductivity. But 
it is also possible that D becomes finite, turning a T = insulator to 
an ideal conductor; for instance, a system of independent particles (e.g. 
one described within a mean field theory scheme), insulating due to the 
band structure, turns to an ideal conductor at T > 0. Finally, it is also 
conceivable that both D and cr^c remain zero at T > 0, a system that 
can be called an "ideal insulator". 

To the above scenaria we should add the possibility that the low fre- 
quency conductivity at finite temperatures is anomalous, e.g. diverging 
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as a power law of the frequency, resulting in an infinite dc conductivity. 
Actually, as we will discuss later (Discussion section), this kind of be- 
havior is fairly common in classical one dimensional nonlinear systems. 

Thus, the first step in characterizing a system is the evaluation of 
the Drude weight at T = in order to find out whether the system 
is conducting or insulating. The peculiarity that has recently been no- 
ticed is that most prototype models, assumed faithful representations 
of the physics of several quasi-one dimensional materials, have finite 
Drude weight also at finite temperatures (even T — > oo), thus imply- 
ing intrinsically ideal conductivity. In other words, interactions do not 
present a sufficient scattering mechanism to turn these systems into nor- 
mal conductors. This behavior is unlike the one observed in the higher 
dimensional version of the same models, that become normal conduc- 
tors at finite temperatures [24] . This unconventional behavior has been 
attributed to the intcgrability of these models. 

To evaluate the Drude weight is not an easy matter as, although 
frequency independent, it represents a transport property and thus it 
cannot be obtained via a thermodynamic derivative (e.g. of the free en- 
ergy). Direct calculation using the optical sum rule eq.(1.6) is obviously 
involved requiring the value of all current matrix elements. A very con- 
venient and physical formulation is the one by W. Kohn, eq.(1.7), that 
generalized at finite temperatures [25] reads, 

D = ^Y.Pn-Q^U^o. (1.10) 

By considering the change of the free energy as a function of flux (that 
vanishes in the thermodynamic limit as it is proportional to the suscep- 
tibility for persistent currents) we can also arrive at an expression for 
the Drude weight as the long time asymptotic value of current-current 
correlations [26], 

D = ^ EPn{n\j\nf = ^{jmt^oo ^ PC,,. (1.11) 

n 

As an example, for a ID tight binding free spinless fermion system with 
nearest neighbor hopping t, the application of a flux (f> modifies the single 
particle dispersion to = —2tcos{k + (f)) giving, 

Do = - sin(7rn) = N{eF)jF 

TT 

D{T)^Do-^{jf (n = i). (1.12) 
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Here, n is the fermion density, N{eF) the density of states and jp the 
current at the Fermi energy. Notice the quadratic decrease with tem- 
perature of the Drude weight that, as we will see later, it is generic even 
for interacting one dimensional fermionic systems out of half-filling. 

In the recent literature, that we will discuss below, the Drude weight 
of integrable systems is evaluated by the BA technique at zero or finite 
temperatures using the Kohn expression (1.10). The difficulty in this 
approach is the need for the estimation of finite size energy corrections 
of the order of 1/L, a rather subtle procedure within this method. 

Another approach, proved particularly efficient in establishing that 
systems with a finite Drude weight at finite temperature exist, uses an 
inequality proposed by Mazur [27] . This inequality states that if a system 
is characterized by conservation laws Qn then: 



Here ( ) denotes a thermodynamic average, the sum is over a subset of 
conserved quantities Qn orthogonal to each other in the sense {QnQm) = 
{Qn)5n,m, = A and we take (A) = 0. 

Thus, for time correlations {A{t)A) with non-singular low frequency 
behavior we can obtain a bound for Caa = ^i^t->oo{A{t)A) , 



For integrable systems, such as the spin-1/2 Heisenberg or Hubbard 
model that are known to possess nontrivial conservation laws because 
of their integrability, useful bounds can be obtained by considering just 
the first non-trivial conservation law. We should stress however that this 
approach has not provided yet a complete picture of the Drude weight 
behavior as we will discuss below in concrete examples. 

Finally, another argument relating the behavior of the Drude weight 
to the (non-) integrability of a model is by the use of Random Matrix 
Theory [28, 29, 25]. It is known that integrable systems are character- 
ized by energy level crossings upon varying a parameter and so Poisson 
statistics in the energy level spacing; thus it can be argued that the 
typical value of diagonal current matrix elements (slope of energy levels 
with respect to an infinitesimal flux) is of the order of one, plausibly 
implying a finite Drude weight according to eq.(l.ll). On the contrary, 
nonintegrable systems, due to level repulsion, are described by Wigner 
or GOE statistics and thus the characteristic value of diagonal current 




(1.13) 




(1.14) 
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matrix elements is of the order of (inversely proportional to the 
density of many body states) implying now a vanishing Drude weight as 
L — ^ oo. 

Besides electrical transport, the thermal conductivity of ID systems 
has recently attracted particular interest; within linear response theory 

it is given by the analogous Grccn-Kubo formula expressed in terms of 
the energy current - energy current dynamic correlation function, 

k{oj) = ^—XjEjE{ijj). (1-15) 

Unlike the conductivity, there is no "mechanical force" (as the flux (j)) 
that can be applied to the system in order to deduce expressions simi- 
lar to the Drude weight, but the long time asymptotic value of energy 
current correlations has an analogous meaning. 

Finally, in magnetic systems, the "spin conductivity" (spin diffusion 
constant) can be probed, for instance, by NMR experiments that mea- 
sure at high temperatures the Fourier transform of spin-spin autocorre- 
lations at the Larmor frequency wjv, 

/ + 00 f 
dt / dqe'^^\Sl{t)St^). (1.16) 

By using the continuity equation, 

u^^{SlSt^)^ = q\jlf_,)^ (1.17) 

for a system where the total spin z— component is conserved, the spin- 
spin dynamic correlations can be analyzed via the corresponding spin- 
current correlations in analogy to electrical transport [30]; the role of 
local charge is played by the z— component of the local magnetization 
(see next section for a more detailed discussion on this point). 

We will now briefly discuss different methods, analytical and numeri- 
cal, that are available for the study of finite temperature dynamic corre- 
lations in strongly interacting systems. Among the analytical approaches 
that have been used for the study of transport and dynamic properties 
of ID systems, each has its own advantages and drawbacks. The tra- 
ditional memory function approach [31] provides a complete picture of 
the temperature/frequency dependence but it is a perturbative method 
based on the assumption of a regular relaxation behavior that might 
be dangerous in ID systems. The high temperature moment expansion 
provides useful information on the possibility of anomalous transport 
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but the extraction of transport coefficients is also based on tfie phe- 
nomenological assumption of regular, diffusive behavior [32]. Progress 
in the exact evaluation of dynamic correlations in integrable systems 
has recently been achieved in the calculation of the Drude weight by 
the Bethe ansatz technique and of the frequency dependent conductiv- 
ity by the form factor method. The Drude weight studies however are 
still controversial as they involve the calculation of finite size corrections, 
while the form factor approach has so far been limited to the calcula- 
tion of zero temperature correlations and mostly in gapped systems. It 
is expected however that progress in BA techniques will provide a full 
picture of the dynamic properties of integrable systems. It is amusing 
to remark the paradoxical situation where the only strongly correlated 
systems for which we can probably have a complete solution of their dy- 
namics are the integrable ones, which however, exactly because of their 
integrability, show unconventional behavior. 

Among numerical simulation techniques, the ED (exact diagonaliza- 
tion) provides exact answers over the full temperature/frequency range 
but of course only on finite size systems [33]. Due to the exponentially 
growing size of Hilbert space, this limits the size of systems that can 
be studied to only about 20 to 30 sites, depending on the complexity 
of the Hamiltonian. We should also remark that, in principle, the full 
excitation spectrum is required for the evaluation of finite temperature 
correlations ^ Furthermore, the obtained frequency spectra are discrete, 
5— functions corresponding to transitions between energy levels, so that 
some ad-hoc smoothing procedure is needed; this is particularly cru- 
cial in attempting to extract the low frequency behavior. Nevertheless, 
finite size scaling in ID systems can provide very useful hints on the 
macroscopic behavior, particularly at high temperatures where all en- 
ergy levels are involved. This regime is the most favorable in attempting 
to simulate the physical situation where the scattering length is less or 
comparable to the system size. 

The Quantum Monte Carlo techniques allow the study of far larger 
systems and they provide directly the dynamic correlations at finite tem- 
peratures but in imaginary time [35]. By analytical continuation, using 
for instance the Maximum Entropy procedure, one is able in principle to 
extract the main features of the frequency dependence; experience shows 
however, that fine issues as the temperature dependence of the Drude 



^In a recent advance, finite temperature dynamic correlations for a prototype model have 
been successfully evaluated using only one quantum state (microcanonical ensemble) [34] . 
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weight or the presence of diffusive behavior which is a low frequency 
property, are difficult to establish reliably. 

Finally, the DMRG method that has been so successful in the study 
of ground state and thermodynamic properties of ID systems, has only 
recently been extended to the reliable study of zero temperature con- 
ductivities in gapped systems [36, 37]. At finite temperatures it is also 
possible to obtain very high accuracy data on autocorrelation functions 
in imaginary time by the use of the transfer matrix DMRG [38]. How- 
ever, similarly to QMC methods, it is very difficult to extract subtle 
information on the finite T dynamics because of the extremely singular 
nature of analytic continuation that hides the useful information even 
for practically exact imaginary time data. 

3. Heisenberg model 

The prototype model for the description of localized magnetism is the 
Heisenberg model. For a one dimensional system the minimal Hamilto- 
nian describing magnetic insulators is, 

H = Y.hi = jj^iSfSUi + SfSl, + A5f 5f+i) (1.18) 
I 1=1 

where Sf (a = x, y, z) are spin operators on site I ranging from the most 
quantum case of spin S=l/2 to classical unit vectors. For S=l/2 the 
system is integrablc by the Bcthc ansatz method and its ground state, 
thermodynamic properties and elementary excitations have well been 
established [10]. As a brief reminder to the discussion that follows, note 
that for J > 0, A > corresponds to an antiferromagnetic coupling while 
A < to a ferromagnetic one; a canonical transformation maps H{A) to 
—H{—A). Further, the anisotropy parameter A describes two regimes, 
the "easy-plane" for |A| < 1 or the "easy-axis" for |A| > 1. The isotropic 
case, occuring in most materials for symmetry reasons, corresponds to 
A = 1. For |A| < 1 the system is gap less and characterized by a linear 
spectrum at low energies, while for A > 1 a gap opens; in particular, 
at A = 1 the elementary excitation spectrum is described by the "des 
Cloiseaux-Pearson" dispersion = -^Isingl. For A < — 1 there is a 
transition to a ferromagnetic ground state. 

In general, other types of terms appear in the description of quasi-lD 
materials STich as longer range or on site anisotropy interactions, but in 
this review we will focus on the prototype model cq.(1.18). 

At this point we should mention that the spin- 1/2 Heisenberg model is 
equivalent to a model of interacting spinless fermions (the "t-V" model) 
obtained by a Jordan- Wigner transformation [39]; 
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L ^11 

H = i-t) + h.c.) + V J2{ni - -)ini+, - -), (1.19) 

1=1 1=1 

where c;(c|^) denote annihilation (creation) operators of spinless fermions 
at site / and n; = c^q. 

The correspondence of parameters is V/t = 2 A and the opening of 
a gap at A > 1 corresponds to an interaction driven metal-insulator 
(Mott-Hubbard type) transition. 

3.1 Currents and dynamic correlations 

Regarding the transport and dynamic properties of the Heisenbcrg 
model, three cases have mostly been discussed: the classical one, the spin 
S=l and the spin S=l/2; the S=l case has been extensively analyzed 
by mapping its low energy physics to a field theory [40] , the nonlinear-a 
model (see section 5). In connection to experiment, the main issue is the 
diffusive vs. ballistic character of spin transport as probed for instance 
by NMR experiments and recently the contribution of magnetic exci- 
tations to the thermal conductivity of quasi-one dimensional materials 

[!]• 

To discuss magnetic transport, we must first define the relevant spin 
and energy currents by the continuity equations of the corresponding 
local spin density Sf (provided the total component is conserved) 
and local energy hf, 

fiQZ 

gives for the spin current, 

/ = E = JT^i^iSli - SfSf^,). (1.21) 
/ I 

Here and thereafter, Va/ = ai — ai-i denotes the discrete gradient of a 
local operator a/. In general (A ^ 0) the spin current does not com- 
mute with the Hamiltonian, [j^,H] ^ 0, so that nontrivial relaxation is 
expected and thus finite spin conductivity at T > 0. 
Similarly, the energy current is obtained by. 
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+ AiSUSfSr^,-Sf_,SfSU,) (1.23) 

We will now briefly comment on the framework for discussing spin 
dynamics and in particular how it is probed by NMR experiments. Ac- 
cording to the spin diffusion phenomenology (for a detailed description 
see ref. [41]) when we consider the {q,uj) correlations of a conserved 
quantity A = J^i^h such as the magnetization or the energy, it is as- 
sumed that it will show a diffusive behavior in the long-time t — oo, 
short wavelength q ^ regime ^. In the language of dynamic correla- 
tion function, diffusive behavior means that the time correlations decay 
as, 

{{Mt),Aom) = 2xaT J ^e'<ii-D^i'\t\ (1.24) 

where Da, XA are the corresponding diffusion constant and static sus- 
ceptibility, respectively. For a ID system, this behavior translates to a 
characteristic dependence of the autocorrelation function. 

Fourier transforming the above expression we obtain, 

SAA{q,i^) = r^dt e-*i({yl,(t),^_,(0)}) ~ ^j^'l^.fl^o - (1-25) 

j-oo ^ \J^Aq ) "T 

By using the continuity equation (1.20), this Lorentzian form can be 
further modified to obtain the current-current correlation function. 



XaDa^^"^ 



SjAjA{q,u) ^ ^ g (1.26) 



which gives the diffusion constant Da by taking the g — limit first 
and then a; — > 0. 

On the other hand, a ballistic behavior is signaled by a 5— function 
form, SjAjA{q,u;) ~ 6{lo — cq), where c is a characteristic velocity of 
the excitations. This (5— function peak moves to zero frequency as g — 
and its weight is proportional to the long time asymptotic of the current- 
current correlations 

CjAjA = SjAjA{q = 0,t^ go). (1.27) 



^This phenomenological statement goes under the name of Ohm's law in the context of 
electrical transport, Fourier's law for heat or Pick's law for diffusion. 
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The above anticommutator correlations are related to the imaginary part 
of the susceptibility x(g,u;), that describes the dissipation, by, 

SAA{q,u) = coi\i{^)^'j,j,{q,u). (1.28) 

In relation to the experimental study of spin dynamics, the NMR has 
developed to a very powerful tool; for instance, the 1/T\ relaxation time 
is directly related to the spin-spin autocorrelation by, 

- ^ / dtcos{uMt){{S!{t\Stm) (1-29) 

where |Ap is the hyperfine coupling [4] and ujn the Larmor frequency. 
Using the relation (1.28), l/Ti gives information (in the high tempera- 
ture limit, Plon — ^ 0) on x"{q,^) as, 

I^TIAI^E^^^- (1-30) 

The diffusive behavior, characterized by the l/\/t decay of the spin cor- 
relations, is extracted in an NMR experiment by analyzing the q ^ 
contribution [5] . It gives a 1/ \/loj^ behavior that is detected as a 1/ y/H 
magnetic field dependence, 

J_ J_ Txjq = 0) 



(1.31) 



considering that the Larmor frequency lon ^ H , Dg being the spin 
diffusion constant and x(q = 0) the static susceptibility. 

3.2 Spin and energy dynamics 

Returning now to the state of spin and energy dynamics, the classical 
Heisenberg model has been extensively studied by numerical simulations, 
the first studies dating from the 70's [42]. Nevertheless, the issue of 
diffusive behavior (even at T = oo where most simulations are carried 
out) still seems not totally clear, the energy and spin showing distinctly 
different dynamics. On the one hand, simulations clearly indicate that 
energy transport is diffusive [43] but on the other hand, the decay of 
spin autocorrelations is probably inconsistent with the expected \/^/t 
law [44, 43] exhibiting long-time tails. 

On the other extreme, for the fully quantum spin S=l/2 model, the 
simplest case is the A = 0, so called XY limit. Here, the spin current 
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commutes with the Hamiltonian resulting in baUistic transport; this can 
also be seen in the fermionic, t-V, version of model that corresponds 
to free spinless fermions (yjt = in eq.(1.19)) where now the charge 
current is conserved. In the infinite temperature limit [fi = 0) the 
spin and energy autocorrelations can be calculated analytically using 
the Jordan- Wigner transformation and are of the form [45]: 

{S!m) = \j^{Jt) (1.32) 

{hi{t)hi) = ^ (j',{Jt) + jUJt)) (1.33) 

which both behave as 1/t for t ^ oo, unlike the l/y/i form in the 
diffusion phenomenology (Jo, Ji are Bessel functions). Actually the /3 = 
limit, often theoretically analyzed for simplicity, is not unrealistic as 
the magnetic exchange energy J can be of the order of a few Kelvin in 
some materials. 



For |A| < 1 the Drude weight at T = has been calculated using the 
BA method [46, 21] and is given by, 

_ TT sin(7r/z/) 

where A = cos^-n/v) ^. For A > 1, D{T = 0) = as the system is 
gapped. 

At finite temperatures, several numerical and analytical studies indi- 
cate that for |A| < 1 the spin transport is ballistic [47, 48, 49, 50, 51], in 
accord with the conjecture that this behavior is related to the integrabil- 
ity of the model [52, 53, 25]. Pursuing this conjecture, one can attempt 
to use the Mazur inequality eq. (1.14) in order to obtain a bound on 
the Drude weight and thus establish that the transport is ballistic. In- 
spection of the known conservation laws for the Heisenberg model [54] 
shows that already the first nontrivial one, Qs, has a physical meaning; 
it corresponds to the energy current, Q3 = and it can be used to 
establish a bound for D [26] , 



/9 {fQz? 
2L (Qi) 



^The parametrization of A in terms of u is common in the BA literature as the formulation 
greatly simplifies for 1/ =integer. 
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This expression can be readily evaluated in the high temperature limit 
{(3^0), 



D{T) > 



P 8A2m2(l/4-m2) 



m = {Sf), 



(1.36) 



2 l + 8A2(l/4 + m2)' 
where m is equal to the magnetization density in the Heisenberg model 
or to n — 1/2 in the equivalent fcrniionic t — V model (n is the density). It 
establishes that ballistic transport is possible at all finite temperatures 
in the Heisenberg {t — V) model; notice however, that the right hand 
side vanishes for m = 0, that corresponds to the specific case of the 
antiferromagnetic regime at zero magnetic field or to the t — V model at 
half-filling. Of course this does not mean that D is indeed zero in these 
cases as this relation provides only a bound. It should also be remarked 
that the obtained bound is proportional to and so we do not recover 
the simple result that D(T) > in the XY-limit. Furthermore, it can be 
shown, using a symmetry argument, that even by taking into account 
all conservation laws the bound remains zero at m = [26]. 

A BA method based calculation of D{T) for |A| < 1 was also per- 
formed [55] , using a procedure proposed for the Hubbard model [56] , that 
relies upon a certain assumption on the flux dependence (see eq.(l.lO)) 
of bound state excitations ("rigid strings"). The resulting behavior is 
summarized in Figs. 1.3 and 1.4. From this analysis the following picture 
emerges: 
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Figure 1.3 -D(A) at var- 
ious temperatures. The 
lowest line is the high 
temperature proportional- 
ity constant Cjj = -D//3. 
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(i) at zero magnetization, in the easy plane antiferromagnetic regime 
(0 < A < 1), the Drude weight decreases at low temperatures as a 
power law D{T) = Dq - const. T°, a = 2/{u - 1); 

(ii) in the ferromagnetic regime, — 1 < A < 0, D{T) decreases quadrati- 
cally with temperature (as in a noninteracting, XY-system); 
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(iii) the same low temperature quadratic behavior is true at any finite 
magnetization; 

(iv) for /3 ^ 0, D{T) = (5Cjj and it can be shown that -^(-A) = D(A) 
by applying a unitary transformation in the expression eq.(l.ll); 

a closed expression for Cjj can be obtained by analytic calculations [57] , 
= (7r/i/-0.5sin(27r/z/))/(167r/z/) for |A| < 1 while Cjj = for A > 1; 

(v) at the isotropic antiferromagnetic point (A = 1), D(T) seems to 
vanish, implying non ballistic transport at all finite temperatures. 



A=+COS(7l/4) 
A=+COS(7l/5) 
A=-COS(7C/4) 

A=-cos(jt/5) 




Figure 1.4 Temperature 
dependence of the Drude 
weight D vs. T [55] 



This last result seems in accord with the most recent NMR data [5]. Of 
course, the low frequency conductivity must also be examined in order 
to determine whether there is no anomalous behavior (e.g. power law 
divergence) that precludes a normal diffusive behavior; such unconven- 
tional behavior is presently debated in classical nonlinear ID systems 
(see final section of Discussion). It should not be surprising if future 
rigorous studies reveal that the isotropic Heisenberg exhibits a singular 
behavior, as it lies at the transition between a gapless and gapped phase. 

In this context, we should also mention that the power law decrease of 
D{T) for < A < 1 is not corroborated by recent QMC simulations [58]. 
The disagreement might be due either to the "rigid string" assumption 
used in the BA analysis or to the very low temperatures, of the order of 
the energy level spacing, that are studied in the QMC simulations ^. 

Considering the limited results obtained so far using the Mazur in- 
equality compared to the exact BA analysis, it remains an open ques- 
tion whether the behavior of the Drude weight can be fully accounted 



^Reliable results for the Drude weight can be obtained by QMC simulations only at low 
temperatures because a sufficiently fine spacing of Matsubara frequencies is required for the 
extrapolation to zero frequency. 
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for solely by a proper consideration of conservation laws present in the 
Heisenberg model. 

For A > 1 numerical simulations [48] and analytical arguments [51] 
indicate that the Drude weight vanishes at all temperatures. In this 
regime, based on ED numerical simulations, it was proposed that a new 
phase might exist, an "ideal insulator", characterized by vanishing Drude 
weight and diffusion constant {dc conductivity in the fermionic version). 
This conjecture remains presently still rather tentative, due to the small 
size of the systems that have been studied so far. 

On the other hand, a semiclassical field theory approach [59] con- 
cluded that gapped systems are diffusive. This approach is based on a 
mapping of the massive excitations to impenetrable classical particles 
of two or more charges (corresponding to different spin directions) that 
propagate diffusively (see section 5) and it has mostly been used for the 
analysis of gapped spin-1 systems. 

In parallel to these developments, the spin S=l/2 Heisenberg model 
was studied in the scaling limit using conformal invariance arguments 
[60, 61]. This field theory approach amounts to considering a linearized 
spectrum and thus neglecting the effects of curvature, a point that we 
will further discuss below in section 5. In particular, it was shown that 
the uniform dynamic susceptibility describes ballistic behavior, the cor- 
responding l/Ti relaxation time was evaluated and the theory was ex- 
tensively compared to experimental data [62]. Notice, however, that a 
later experimental NMR work [5] concludes that the q = mode of spin 
transport is ballistic at the T = limit, but has a diffusion-like contri- 
bution at finite temperatures even for T << J. We should remark that, 
over the years, the most common interpretation of NMR experiments 
was within the diffusion phenomenology, as for instance for the S = 5/2 
TMMC compound [63]. 

Finally, the finite {q,oo) response functions of the S=l/2 model at 
r = were studied by the bosonization technique [64] after mapping it 
to spinless fermions (eq.(1.19)). For A < 1, the conductivity shows the 
typical ballistic form; for A > 1 it vanishes below the gap, showing a 
square-root frequency dependence above. 

Turning now to energy transport, it is easy to see that the energy 
current is a conserved quantity [65, 26] for all values of the anisotropy A 
implying that the currents do not decay and so the thermal conductivity 
is infinite. This peculiarity has also been noticed by an earlier analysis 
of moments at infinite temperature [66]. So the quantity characterizing 
thermal transport is the equal time correlation (j^j^) that represents 



20 



> 
■a 

o 

O 

o 

'rt 2 



' 1 


' 1 ' 


















7^/3 




/ /' ^ 
1 /' ^ 




— Y=Ji/4 


- 




\ 


— Y=7r/5 




/ ■/ / 


■ \ 

■«■ N ^. \ 


y=TZl6 




//■// 


^ ■ \ 

^' ^\ 


— y=Q 




ilf 
/// ' -■■ 
/'// ,•■ 


N N ■ \ 






//// ■■ 
/,// .■ 
/,/' -■■ 








I'l ■■ 

— Ill :' 

hi ; 








i; : 
















.' 








ll: / 








S: ' 
ft' 

f , 1 


1 







0.5 



1 

T/J 



1.5 



Figure 1.5. Thermal conductivity, {j^ j^) in units of J^, for various anisotropy pa- 
rameters A = cos (7) [67]. 



the weight under the low frequency peak that will develop from the zero 
frequency (5— function when a dissipative mechanism is introduced. This 
picture is analogous to that of the electrical conductivity illustrated in 
Fig. 1.1. It implies that, given an estimate of the temperature depen- 
dence of the characteristic scattering time one is able to extract the 
value of the dc thermal conductivity, further assuming some form (e.g. 
eq.(1.8)) for the low frequency behavior. 

This quantity has also recently been exactly calculated using the BA 
method [67]; it is shown in Fig. 1.5. 

Furthermore, the experimental observation of unusually high thermal 
conductivity in ladder compounds [1] motivated the theoretical study of 
the thermal Drude weight in ID anisotropic, frustrated and ladder spin- 
1/2 systems [68, 69]; the proposal of unconventional thermal transport 
in these systems is still debated. 

Finally, the S=l Heisenberg chain shows a qualitatively different be- 
havior characterized by the presence of an energy (Haldanc) gap at low 
energies. The S=l Heisenberg model is not integrable but the physics 
at low energies is usually mapped onto the quantum nonlinear-cr model 
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that is again an integrable system. The results known on this model will 
be briefly discussed in section 5 along with a semiclassical approach to 
describe this type of gapped systems. The same low energy mapping is 
used for the analysis of "ladder" compounds. 

As a guide to experimental investigations and theoretical studies, we 
can recapitulate the above discussion of the dynamics of the Heisenberg 
S=l/2 model as follows. It seems clear that ballistic behavior at all 
temperatures should be expected in the easy-plane regime and at all fi- 
nite magnetizations, while the isotropic point is a subtle borderline case. 
The behavior in the easy-axis antiferromagnetic regime might be par- 
ticularly interesting and it is not settled at the moment. Exceptionally 
high thermal conductivity should be expected in all regimes. 

To complete the above picture, we should stress that not much is 
known on the low frequency behavior of the conductivities at finite tem- 
perature. This leaves open the possibility of unconventional behavior, 
neither ballistic nor simple diffusive but one characterized by long time 
tails, giving rise to power law (or logarithmic) behavior at low frequen- 
cies. 

4. Hubbard model 

The prototype model for the description of electron-electron correla- 
tions is the Hubbard model given by the Hamiltonian, 

H = Y.hi = i-t) + h.c.) +UJ2 ni^m (1-37) 

/ (7,1 I 

where cicr{c]^) are annihilation (creation) operators of fermions with spin 
=1: i site I and n^o- = c^^ci^- 

At half-filling (n=l, 1 fermion per site) it describes a Mott-Hubbard 
insulator for any value of the repulsive interaction U > 0, while it is a 
metal at any other filling. 

The one dimensional Hubbard model is also integrable by the Bethe 
ansatz method and its phase diagram, elementary excitations, correla- 
tion functions have been extensively studied [10, 70]. 

4.1 Currents 

Similarly to the Heisenberg model, we can discuss the electrical, spin 
and thermal conductivity by defining the charge j, spin and energy 

currents from the respective continuity equations of the local particle 
density ni, 
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^ + Vji =0, j = Y^ji = Y^ji, = i-t) ^(i4q+i, + h.c), (1.38) 
I la a,l 

spin density n;| — n^, 

and energy density hi, 

f^ + Vjf = 0, f = Y.jFa (1-40) 

jFa = i-t)'^ii4^ia'^i-ia + h.c.)-^ji,„{ni-„ + ni+i-^-l). 

4.2 Electrical and thermal transport 

With respect to the electrical conductivity the interaction U and den- 
sity dependence of the Drude weight D at zero temperature has been 
established using the BA method [71, 72, 73] (see Fig. 1.6). There are 
two simple limits: 

(i) The free fermion case [7 = where j is conserved and Dq = ^ sin ^ 
where n is the density of fermions {n = 2kp/n). Here Dq vanishes for 
an empty band n = and a filled band n = 2, being maximum at half 
filling, n = 1. 

(ii) Another simple limit is ?7 = cxd. Since in this case the double occupa- 
tion of sites is forbidden, fermions behave effectively as spinless fermions 
and the result is I?o = |:| sin(7rn)|; here Dq vanishes also at half filling. 

Analytical results in ID indicate that the Dq = value at half fill- 
ing persists in the Hubbard model for all U > 0, whereby the density 
dependence D{n) is between the limits U = and U = oo. The insulat- 
ing state at half filling is a generic feature of a wider class of ID models 
characterized by repulsive interactions, such as the t-V model (discussed 
above), the t-J model etc. 

In Fig. 1.6, along with the Drude weight, the zero temperature (ballis- 
tic) Hall constant Rh of a quasi-lD system is also shown. According to 
a recent formulation [74], Rh can be expressed in terms of the derivative 
of the Drude weight with respect to the density. 



1 dD 
~ ~D~dn' 



(1.41) 




Figure 1.6. Drude weight D and Rh for the quasi-lD Hubbard model from expression 
(1.41). 

The Hall constant is the classical way for determining the sign of the 
charge carriers. For a strictly one dimensional system of course it makes 
no sense to discuss the Hall effect; but if we consider a quasi-one di- 
mensional system with interchain coupling characterized by a hopping 
t' 0, then within this formulation we recover a simple picture for the 
behavior of the sign of carriers as a function of interaction. In agree- 
ment with intuitive semiclassical arguments, the Hall constant behaves 
as Rh — —1/n at low densities changing to Rh — +1/5(5 = 1 — n) near 
half-filling, with the turning point depending on the strength of the in- 
teraction U. Notice that if D oc n with a small proportionality constant, 
that would be interpreted within a single particle picture as indicative 
of a large effective mass, then we would still find Rh — —1/n. This 
observation might be relevant in the context of recent optical and Hall 
experiments [75, 7] where a small Drude weight is observed although the 
Hall constant indicates a carrier density of order of one. 

Recently, using the form factor and DMRG methods the frequency 
dependence of the conductivity at half filling and at T = has also been 




Figure 1.7. Optical conductivity at T = for U/t — 3,6,12 (from left to right) 
calculated with DDMRG on a 128-site lattice [76]. Inset: a{oj) for U/t = 12 (dashed) 
and 40 (solid) calculated on a 64-site chain. 



studied [76] and is shown in Fig. 1.7. The DMRG method provided 
the entire absorption spectrum for all but very small couplings where 
the field theoretical approach was used; the two methods are in excel- 
lent agreement in their common regime of applicability. As expected, 
the Drude weight is zero, signaling an insulating state (for a detailed 
analysis of the scaling of D with system size at and close to half-filling, 
see [77]) and the finite frequency conductivity vanishes up to the gap. 
Above the gap, a square root dependence is observed btit not a diver- 
gence; this behavior is in contrast to that obtained by the Luttinger 
liquid method [78] and it is typical of a Peierls (band) insulator where 
a divergence occurs. This absence of a singularity is also in agreement 
with a rigorous analysis of the sine-Gordon (sG) field theory (see sec- 
tion 5), the generic low energy effective model for the description of a 
Mott-Hubbard insulator. 

To complete the zero temperature picture, the frequency dependent 
conductivity of the Hubbard model out of half-filling has been studied 
using results from the BA method and symmetries [79]. A broad ab- 
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sorption band was found separated from the Drude peak at w = by a 
pseudogap; this pseudo-gap behavior is in contrast to the o;^ dependence 
found within the Luttinger hquid analysis [78]. 

Again, at all finite temperatures the transport is ballistic characterized 
by a finite Drude weight. In an identical way to the Heisenberg model, 
this can easily be established by the Mazur inequality using the first 
nontrivial conservation law Qs- For the Hubbard model differs from 
the energy current by the replacement of U by U/2 [26]. Evaluating 
{jQzfKQl) for ^ ^ we obtain, 



Dm > JL^J9^ = /3 [L^E.2n.(l-n.)(2n-.-l)]^ 
^ '-2L {QD 2Ea2n,(l-n<,)[l + C/2(2n2_,-2n_, + l)]' 

(1.42) 

where ricr are the densities of a =t, J, fermions. 

By inspection we can again see that from this inequality we cannot 
obtain a finite bound for D{T) for n| + nj = 1. Nevertheless, a full 
BA calculation [56] seems to show that the Drude weight at half-filling 
is exponentially activated D{T) ~ --^e~^f"p/^ at low temperatures and 

decreases as out of half-filling. Thus the zero temperature insulator 
turns to an ideal conductor at finite temperatures. Notice that this be- 
havior is different from the one in the Heisenberg (or "t-V") model in 
the gapped phase (A > 1) where the Drude weight seems to vanish at 
all finite temperatures. We can conjecture that this distinct behavior 
of insulating phases can be understood in the framework of the corre- 
sponding low energy sine-Gordon field theory as these two models map 
to different parameter regimes of the sG model [64, 13]. A very similar 
calculation, using the Mazur inequality, can also be carried out for the 
long time asymptotics of the spin current, f , correlations. It gives a 
finite bound, and thus ballistic spin transport for n-^ — n[ ^ 0; no BA 
calculation has so far been performed for the spin conductivity. 

On the thermal conductivity we find similar results, namely a finite 
value on the long time decay of energy current correlations, which can 
readily be evaluated for /3 — > [26], 



.li^OW) = C«.>^. (1.43) 

Again this inequality gives a finite bound for a system out of half- filling 
as long asn^+ni 7^ and this for any magnetization. For this model the 
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actual temperature dependence of Cjeje = limj^oo(j^(i)j^) = Cjeje 
has not yet been evaluated. Finally, the low temperature thermoelectric 
power was studied using the Bethe ansatz picture for the charge (holons) 
and spin (spinous) excitations [80]. The resulting sign of the ther- 
mopower close to the Mott-Hubbard insulating phase is consistent with 
the one derived from the Hall constant above, S ~ sign{l — n)T\m*\/\l — 

In summary, we have shown that the prototype model for describing 
electron correlations in one dimensional systems, the Hubbard model, 
shows unconventional, ballistic charge, spin and thermal transport at 
all finite temperatures. Of course real quasi-one dimensional materials 
are presumably characterized by longer range than the Hubbard U in- 
teractions. So, although the above picture should be taken into account 
in the interpretation of experiments, (quasi-) one dimensional magnetic 
compounds might presently appear as better candidates for the exper- 
imental observation of these effects. Theoretically, the full frequency 
dependence of the conductivities at finite temperatures remains to be 
established. 

5. Effective field theories 

An alternative to analyzing the transport of quasi-one dimensional 
materials within microscopic models, as described in previous sections, 
is to approach the problem within effective low energy models for in- 
teracting electrons, i.e. starting with the Luttinger liquid Hamiltonian. 
This path is very attractive since it represents the counterpart of the 
usual Landau phenomenological approach to Fermi liquid in higher-D 
electronic systems. It should be pointed out that even in a 3D system 
the continuum field theory is not enough to describe a current decay and 
Umklapp processes are finally responsible for a finite intrinsic resistivity 
p(r) ocT^ [81]. 

In an effective (low energy) field theoretical model for ID interacting 
electrons the band dispersion around the Fermi momenta k = ikp is 
linearized and left- and right- moving excitations are defined. Apart 
from Umklapp terms, the model of interacting fermions can then be 
mapped onto the well known Luttinger liquid Hamiltonian [16, 13] and 
analyzed via the bosonization representation. In particular one obtains 
for the charge sector, 

Ho = ^J dx[upKp{7TU,f + ^(9,</.p)'], (1.44) 

where the charge density is p(x) = dx(pp and Hp is the conjugate momen- 
tum to (f)p. Interactions appear only via the velocity parameter Up and 
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Luttinger exponent Kp. The charge current in such a Luttinger model, 
j = \/2upKpIlp, is clearly conserved in the absence of additional terms. 
Umklapp terms can as well be represented with boson operators, 

Hj_ = gj_ / dxcos{mV8(pp{x) + 6x), (1-45) 

where m is the commensurability parameter (m = 1 at half-filling - one 
particle per site, m = 2 for quarter filling - one particle for two sites etc) 
and S the doping deviation from the commensurate filling. In principle, 
the mapping of a particular (tight binding) microscopic model onto a 
field theory model, e.g. via perturbation theory, generates terms H„i 
with arbitrary m. While Umklapp terms arc irrelevant in the sense of 
universal scaling of the static properties, they appear to be crucial for 
transport. They drive a metal at half-filling to an insulator, while at an 
arbitrary (incommensurate) filling they should cause a finite resistivity 
since the current is not conserved any more (for an overview of the 
transport properties emerging within the Luttinger liquid picture see 
[82]). 

However, the proper treatment of transport within the Luttinger pic- 
ture in the presence of Umklapp processes is quite involved and even 
controversial. Giamarchi [78] first calculated the cff"cct of Umklapp scat- 
tering within lowest order perturbation theory for the memory function 
M(u;); he thus determined the low-w behavior of the dynamical conduc- 
tivity a{uj) oc l/{uj + M{lo)) that yielded a nonzero finite temperature 
conductivity. At the same time he realized, by using the Luther-Emery 
method, that the Umklapp term can be absorbed in the Hamiltonian 
in such a way as to conserve the current and pointed out the possi- 
bility of infinite dc conductivity even in the presence of Umklapp. A 
similar lowest-order analysis [83] for general commensurate filling pre- 
dicts at T = that (7{io) oc uj'^~'^ and the resistivity p{T) oc T'^ with 
u = An^Kp — 3. On the other hand, Rosch and Andrei [84] pointed 
out that even in the presence of general Umklapp terms there exist par- 
ticular operators, linear combinations of the translation operator and 
number difference between left- and right- moving electrons, which are 
conserved. Since in general such operators have a nonvanishing over- 
lap with the current operator j, this leads to finite Z?(T > 0) > if 
only one Umklapp term is considered. At least the interplay of two 
noncommuting Umklapp processes is needed to yield a finite resistivity 
p(T > 0) > 0. 

From a different perspective Ogata and Anderson [85] argued that 
because of spin-charge separation in ID systems an effect analogous to 
phonon drag (in this case spinon-holon drag) appears that leads to a fi- 
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nite dissipation. Using a Landauer like semi-phenomenological approach 
they conchidcd the existence of a hnear-T resistivity and hnear frequency 
dependence of the optical conductivity. 

The bosonization of the Luttinger liquid model leads [13] to the quan- 
tum sine-Gordon model (eq.(1.45)) which is an integrable system and has 
extensively been studied as a prototype nonlinear quantum (or classi- 
cal) field theory. It is the generic field theory for describing the low en- 
ergy properties of one dimensional Mott insulators. The thermodynamic 
properties and excitation spectrum consisting of solitons/antisolitons 
and breather states have been established by semiclassical and BA tech- 
niques [10] . Presently, there is an effort to determine the transport prop- 
erties of this model rigorously. In particular, the frequency dependence 
of the zero temperature conductivity in the commensurate (insulating) 
phase, zero soliton sector, has been evaluated using the form factor ap- 
proach [86]. The main result is that the square root singularity at the 
optical gap, characteristic of band insulators, is generally absent and 
appears only at the Luther-Emery point; furthermore, the perturbative 
result [78] is recovered only at relatively high frequencies. Besides these 
studies, the Drude weight and optical response near the metal-insulator 
transition, in the incommensurate phase at zero temperature, have also 
been studied by Bethe ansatz [87] and semiclassical methods [88]. Still, a 
rigorous evahiation of the Drude weight and frequency dependence of the 
conductivity at finite temperatures is missing; nevertheless, we can plau- 
sibly argue that because of the integrability of the sine-Gordon model, 
it will turn out that also this model describes an ideal conductor at 
least over some interaction range. Thus, it might remain an open ques- 
tion which scattering processes and/or band curvature must be taken 
into account in order to recover a normal, diffusive behavior at finite 
temperatures. 

Finally, it is well known [89] that the spectrum of integer spin and 
even-leg ladder systems is gapped and that the low energy physics is 
described by the one-dimensional quantum 0(3) nonlinear sigma model 



In imaginary time r the action at inverse temperature (3 is given by 



where x is the spatial coordinate, c a characteristic velocity, a = 1,2,3 
is an 0(3) vector index and na{x,T) a unit vector field n^(a;,r) = 1. 

In a series of works, Sachdev and collaborators [90, 91, 92] developed a 
picture of the low and intermediate temperature spin dynamics based on 



[40]. 




(1.46) 
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the idea that the spin excitations can be mapped to an integrable model 
describing a classical gas of impenetrable particles (of a certain number 
of species depending on the spin), a problem that can be treated ana- 
lytically. Within this framework they have extensively analyzed NMR 
experiments on S=l compounds [93] and they concluded that these sys- 
tems behave diffusively. In contrast to this semiclassical approach, using 
the Bethe ansatz solution of the quantum nonlinear— a model [94], Fu- 
jimoto [95] found a finite Drude weight, exponentially activated with 
temperature, and he thus concluded that the spin transport at finite 
temperatures is ballistic. The origin of this discrepancy is not clear at 
the moment and can be due either to a subtle role of quantum effects 
on the dynamics that is neglected in the semiclassical approach or to a 
particular limiting procedure (the magnetic field going to zero) in the 
BA solution. 

6. Discussion 

We hope that the above presentation demonstrated that the transport 
theory of one dimensional quantum systems is a rapidly progressing field, 
fueled by both theoretical and experimental developments. Still, on the 
question, what is the finite temperature conductivity of bulk electronic 
or magnetic systems described by strongly interacting one dimensional 
Hamiltonians, it is fair to say that no definite answer has so far emerged 
nor there is a clear picture of the relevant scattering mechanisms. 

In this context, it is interesting and instructive to draw an analogy 
with the development of the respective field in classical physics, namely 
the finite temperature transport in one dimensional nonlinear systems. 
Interestingly, in this domain we are also witnessing a flurry of activity 
after several decades of studies. Again, the issue of ballistic versus diffu- 
sive (usually energy) transport in a variety of models and the necessary 
ingredients for observing normal behavior is sharply debated [96, 97]. 
Similarly to the quantum systems, numerical simulations are intensely 
employed along with analytical approaches and discussions on the con- 
ceptual foundations of transport theory. 

For quantum systems it is reasonable to expect that the finite temper- 
ature transport properties of integrable models will, in the near future, 
be amenable to rigorous analysis by mathematical techniques, for in- 
stance in the framework of the Bethe ansatz method. At the same time, 
as we mentioned earlier, it is amusing to notice that the integrable sys- 
tems that we can exactly analyze, present singular transport properties 
presumably exactly because of their integrability. 
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To obtain normal behavior, it is reasonable to invoke perturbations 
destroying the integrability of the model, as for instance longer range in- 
teractions, interchain coupling, coupling to phonons, disorder etc. In this 
scenario, it is then necessary to find ways to study the eflcect of pertur- 
bations around an integrable system and in particular to determine the 
vicinity in parameter space around the singular-integrable point where 
unconventional transport can be detected. This issue is also extensively 
studied in classical systems as it is the most relevant in the interpre- 
tation of experiments and in estimating the prospects for technological 
realizations. It is worth keeping in mind the possibility that integrable 
interactions actually render the system more immune to perturbations, 
an effect well known and exploited in classical nonlinear systems [9]. 

Related to this line of argument is the following question. If integrable 
models show ballistic transport and low energy effective theories like 
the sine-Gordon model are also integrable, then which mechanisms are 
necessary to obtain dissipative behavior ? 

Of course it is also possible that the conventional picture according 
to which only integrable systems show ballistic transport might well be 
challenged. One dimensional nonintegrable quantum systems could also 
show singular transport in the form cither of a finite Drude weight or 
low frequency anomalies. This behavior has been observed in classical 
nonintegrable nonlinear systems where the current correlations decay 
to zero in the long time limit but too slowly, so that the integral over 
time (giving the dc conductivity) diverges. The opposite behavior might 
also be realized, namely that integrable quantum systems show normal 
diffusive transport in some region of interaction parameter space (this 
possibility was raised in the case of gapped systems as the easy-axis spin 
1/2 Heisenberg model or the quantum nonlinear— o" model, see section 
5). Furthermore, the issue of the crossover of the dynamics between 
quantum and classical systems has, at the moment, very little been 
explored and in particular the question whether quantum fluctuations 
might stabilize ballistic transport behavior. 

To address all the above open issues there is a clear need for the 
development of reliable analytical and numerical simulation techniques 
(as the DMRG or QMC) to tackle the evaluation of dynamic correla- 
tions at low temperatures. In particular, progress is needed to include 
the coupling between the diflFerent, magnetic, electronic and phononic, 
excitations. 

In summary, one of the most fascinating aspects in this field is to 
understand the extent to which the so successful physics, experimental 
and technological realizations of classical (integrable) nonlinear systems 
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can be carried over to the quantum world of many body (quasi-) one 
dimensional electronic or magnetic strongly interacting systems. This 
effort is accompanied by the experimental challenge to synthesize novel 
materials/systems that realize this physics. 
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